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A  DIFFERENCE  SCHEME  WITH  AUTOCONTROL  ARTIHCIAL  VISCOSITY 
TO  PREDICT  ABLATED  NOSETIP  SHAPE 


Yang  Maozhao 

(China  Aerodynamic  Research  and  Development  Centre) 


Abstract;  The  trio-layer  explicit  difference  scheme  which  is  added 
artificial  viscosity  item  is  found  by  using  filteration  in  the 
computation  of  ablated  nosetip  shape.  The  scheme  is  of  first  order 
accuracy  in  regions  where  shape  change  is  wavy  and  second  order 
accuracy  in  regions  where  is  Smooth.  Numerical  experiments  show 


that  the  scheme  developed  is  effective.  ^ 

Subject  Terms:  S\.blative  nose  cone,  Reent/^^eatifrg,  Aerodyjjamic' 
shape,  Finite  difference  theoryv 
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n  Introduction  /]; 


During  reentry  the  aerodynamic  heating  is  very  strong  in 
region  of  missile  nosetip  where  the  ablative  rate  of  surface  area  is 


very  large  causing  great  shape  change.  The  prediction  of  ablated 
nosetip  shape  has  thus  been  a  crucial  problem  in  ablative  coupling 
calculation.  Some  people  have  used  eigen  linear  methods  to  solve 


shape  equations  and  obtained  the  dispersion  data  of  shape  versus 
time  The  data  have  then  been  processed  by  three  sampling 


smoothing  and  finally  given  the  tip  shape  at  each  calculation  time. 
Other  people  have  solved  shape  equations  by  immediate  application 
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of  difference  theory  ^  \  To  remove  the  data  fluctuation  in 
calculations  .they  added  artificial  viscosity  item  to  difference  schemes 
appropriately.  These  methods  have  been  successful  in  region  where 
shape  change  is  not  great.  However,  in  case  some  ablated  shapes 
such  as  deep  hollow  and  long  nose  cone  happen,  these  calculations 
may  become  difficult  or  may  lack  accuracy  at  some  specific 
points.  In  order  to  solve  this  problem,  we,  inspired  by  the  shock 
wave  data  processing  have  tried  auto  Alteration  of  shape 

parameters  in  calculation  of  ablated  reentry  nosetip  shape,  found  a 
trio-layer  explicit  difference  scheme,  solved  shape  equations  by 
difference  theory  and  successfully  calculated  the  types  of  shapes 
with  deep  hollow  and  long  nose  cone.  Numerical  experiments  have 
showed  that  this  scheme  is  effective  in  guaranting  required 
calculation  accuracy  as  well  as  enlarging  time  steps .  It  has  also  been 
proved  that  this  scheme  is  of  first  order  accuracy  in  regions  where 
shape  change  is  wavy  and  it  is  of  second  order  accuracy  in  region 
where  shape  change  is  Smooth. 

2.  Difference  Scheme  with  Autocontrol  Artificial  Viscosity  for  Shape 

Prediction 

Figure  1  shows  the  coordinate  system  (x,  r)  to  which  a  missile 

tied  where  the  two  curves  represent  nosetip  shapes  at  t;rne  t  and 

j-(x,t) 

t-i-At,  respectively,  and  the  surface  inclination  is  tg  6p  =  — —  .  Let 

V.oo  be  the  ablative  velocity  along  the  surface,  it  is  easy  to  derive 

the  deforming  velocity  of  shape  along  x  direction  from  Figure  1 
which  is 


r 


KEY:  l--shape  at  time  t;  2--  shape  at  time  t  +  At 

Figure  1  Coordinate  system  for  shape  calculation 
—  sinp. 

or  in  other  words  the  differential  equation  for  ablated  nosetip  shape 
can  be  written  in  the  form  of 


Taking  the  derivative  of  (1)  with  respect  to  r  gives  the  differential 

dx 

equation  for  velocity  of  surface  inclination  change.  Let  ^ 


we  get 


=  ViV(v 

^  l+y^3r  ar  ’ 


(1)  and  (2)  constitute  a  set  of  nonlinear  partial  differential  equations 
which  can  be  solved  by  developing  some  difference  scheme.  After 
that  shape  parameters  are  smoothed  by  three  sampling  techniques. 
Our  paper  introduces  a  trio-layer  explicit  difference  scheme  with 
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autocontrol  artificial  viscosity  item.  Combining  (1)  and  (2)  as  a  set, 
we  solved  it  by  difference  theory  without  the  need  of  three  sampling 
evening  up.  For  that,  first  partially  linearize  equation  (2).  Let 

a  = 

A  =  V.«,y/C 

Equations  (1)  and  (2)  then  become 

0X 

ir=“  <3) 

3y  dy 

—  =A— +p  (4) 

Assuming  that  a,  p  and  A  are  temporarily  fixed,  (3)  and  (4)  then 
become  linear  equations.  Following  paper  [4]  to  induce  filtration 
function  Q,  the  equations  can  be  solved  by  developing  trio-layer 
explicit  difference  schemes.  That  is, 

Xj  -Kl-Qi)Xj  +— (Xj+i  +  Xj.,  )1 


"  +  '  ,  n-l  ^2  ,  n-I  n-1 

yj  -l(i-Q2)yj  +-2-<yj+i  +  yj.i  )l 


dy  n  n 
=  (A^)j  .  p.  (6) 
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or  in  other  rewritten  form, 
n  +  1  n - 1 


n  Ql  Ar" 

'  =  a-  - 

J  4  At 


n- 1  n- 1  n-1 

+>^-1 


n+1  n-1 


,,^,n  n  Q2 
'^ar  4  At' 


n-1  n-1  n-1 

yj+,-2yj  +yj., 


It  is  obvious  that  Ar->0  at  At^O,  so  the  differential  equations 
obtained  by  revising  equations  (7)  and  (8)  are 

ax  Ql  Ar^  a^x 

—  =a  +  -r- - TT  (9) 

at  4  At  ar2 


ay  ay  Ql  Ar^  a^y 


(10) 


Equations  (9)  and  (10)  have  the  additional  artificial  viscosity  item 

a^x  a^y 

V] — 2*  ^nd  V2 — ~  compared  to  their  old  equations  (3)  and  (4)  where 

ar"^  ar*" 

Q  Ar^  Q2  Ar^ 

V|=~-“^^and  V2=  “•  .  Their  magnitude  is  closely  related  to 

filteration  function  Q.  vare  second-order  infinitesimals  in  case  that  Q 
is  a  first-order  infinitesimal,  that  is,  Q^0(At,Ar),  they  are  first-order 

infinitesimals  in  case  that  Q  is  a  zeroth-order  infinitesimal,  that  is, 
Q^O(l).  The  filt  ration  function  Q  is  defined  as  a  shape  parameter  x, 

the  relative  fluctuation  value  of  y,  that  is 


(11) 


where  Kj  and  K2  are  constants  usually  taking  values  between  0  and 
1.  Due  to  the  definition  of  Q,  in  regions  where  shape  change  is  wavy 

the  numerators  and  denominators  in  equations  (11)  and  (12)  may 
reach  the  same  numerical  order  at  Q^O(l),  so  the  fore-mentioned 

difference  scheme  is  of  first  order  accuracy.  In  regions  where  shape 
change  is  Smooth,  the  numerators  in  equations  (11)  and  (12)  are 

second-order  infinitesimals  while  the  denominators  are  zeroth-order 
or  first-order  infinitesimals,  so  we  have  at  least  Q...,0  (At,Ar), 

therefore  the  fore-mentioned  difference  scheme  is  of  second  order 


accuracy.  It  is  not  difficult  to  prove  that  this  scheme  is  stable 

provided  the  Alteration  function  in  difference  equations  (7)  and  (8) 
^  2  1 

Qe  [0,^  1  ^  large  number  of  numerical  experiments  show  that 


calculation  time  stens  can  be  enlarged  provided  the  required  accuracy 
is  guaranteed.  Being  compared  with  previous  calculations  the  time 
can  be  reduced  by  one  third  in  computersalong  a  reentry  trajectory. 

During  detailed  calculation  we  have  studied  the  cases  of 
equistep  network,  changeable  step  network  and  changeable  time 
step  along  x  and  r  directions,  respectively.  Besides,  for  comparison 
we  delveloped  several  schemes  including  second  difference  scheme 
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smoothed  by  sampling  and  first  difference  scheme  smoothed  by 
sampling  which  are  discussed  briefly  in  the  next  section. 


3.  Calculation  Results  and  Comparison  Analysis 
Shape  change  of  ablated  reentry  nosetip  is  very  significant 
according  to  surface  roughness  of  missile  nosetip  and  mechanical 
erosion  of  material.  Using  difference  schemes  with  autocontrol 

artificial  viscosity  item  to  solve  ablative  shape  equations  using 
difference  along  reentry  ballistic,  we  have  successfully  calculated  the 
shape  upon  landing  with  deep-hollows  and  long  nose  cone.  Figure  2 
shows  the  time  history  of  ablative  shape  change  obtained  by  our 
scheme  calculation.  Due  to  this  figure  the  final  shape  data  after 
autocontrol  filteration  scheme  processing  remain  smooth  and 

continuous  even  in  the  upper  region  where  shape  shange  gradient  is 
significant.  It  makes  ablative  shape  coupling  calculation  more  simple 
and  concise  without  necessity  of  sampling  smoothing  and  evening  up. 
Figure  3  gives  shape  upon  landing  calculated  by  different  difference 
schemes.  Comparing  the  curves  in  the  figure,  it  can  be  found  that  the 
scheme  with  artificial  viscosity  is  quite  similar  to  the  shape 
calculated  by  first  order  scheme  at  the  upper  region  because  both 
these  schemes  are  of  first  order  accuracy  in  this  region.  The 
calculated  shapes  differ  significantly  from  each  other  in  the  shoulder 
region  where  the  two  schemes  are  different  in  accuracy.  Shape  (2)  is 
rather  close  to  shape  (3)  becuase  they  are  both  of  second  order 
accuracy.  Three  schemes  yield  basically  the  similar  sharp  nose  cone 
shape  at  upper  region  which  means  that  at  upper  region  surface 
roughness  has  a  leading  effect  in  controlling  shape  change  while  the 
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effect  of  diversity  of  difference  schemes  on  shape  calculation  is  less 


Figure  2  The  history  ^of  shape ^(^*ange  of  ablated  nosetip 


Figure  3  The  shape  upon  landing  calculated  using  different 

difference  schemes 

1 — first  order  scheme  and  sampling;  2 — second  order 
scheme  and  sampling;  3 — artificial  viscosity  scheme 
4 — initial  shape 


KEY  : 


Figure  4  shows  effect  on  shape  calculation  due  to  different 
directions  in  network  calculations.  In  artificial  viscosity  scheme,  axial 
equistep  network  and  radial  equistep  network  have  little  difference 
in  effect  on  shape  calculation,  however,  radial  equistep  network  can 
make  the  calculated  points  crowding  together  near  stationary  point 
region  provided  the  points  before  tip  shoulder  point  have  the  same 
count.  Therefore,  when  the  inflection  point  moves  towards 
stationary  point,  prediction  becomes  more  accurate,  resulting  in  more 
meticulous  calculation  of  pressure  and  thermal  flow  near  stationary 
point 


Figure  4  Effects  on  shape  calculation  due  to  different  network 

methods 

KEY:  4--initial  shape;  l--x  equistep  artificial  viscosity 

scheme;  2--r  equistep  artificial  viscosity  scheme; 
3--X  equi  second  order  scheme  and  sampling. 
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